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Abstract 

A variational u = 2/2> state, which unifies the sharp edge picture of Mac- 
Donald with the soft edge picture of Chang and of Beenakker is presented 
and studied in detail. Using an exact relation between correlation functions 
of this state and those of the Laughlin v = l/'i wavefunction, the correlation 
functions of the z/ = 2/3 state are determined via a classical Monte Carlo 
calculation, for systems up to 50 electrons. It is found that as a function of 
the slope of the confining potential there is a sharp transition of the ground 
state from one description to the other. This transition should be observable 

in tunneling experiments through quantum dots. 
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I. INTRODUCTION 



The possibility of observing the fractional charge of edge states in the fractional quantum 
Hall regime has generated considerable excitement in recent years |]1| . While some indications 
of the fractional charge of the edge have been reported , the situation is far from resolved. 
In contrast, in the integer regime there is a profusion of information on the dependence of 
the edge states on the number of electrons, the electrostatic potential (gate voltage) and the 
magnetic field, obtained mainly from tunneling measurements through quantum dots 0. 
This motivated theoretical investigations into the question of tunneling into a quantum dot 
in the fractional quantum Hall regime [P,^. 

Recently it was argued that tunneling into a g = ±e/m edge state will reduce the 
tunneling amplitude by a factor of A^~(™-~i)/2 relative to the integer case. Hence, tunneling 
measurements through a small system in the fractional quantum Hall regime indeed offer 
the possibility of directly probing the composition of the edge structure of the system. At 
zero or low magnetic fields the conductance consists of a series of well separated peaks 
0, each corresponding to an additional electron added to the system. If the average peak 
conductance can be studied as a function of the number of electrons in the system, the 
theory predicts that for filling factor v = 1/3, for example, the peak amplitudes will fall 
as Since, however, tunneling peaks are observed when there is already a substantial 

number of electrons in the system, the tunneling amplitude in this regime may already fall 
below experimental sensitivity. 

This motivated the suggestion that it would be advantageous to study the z/ = 2/3 
regime, where the edge structure is particularly intriguing, and several theories have been 
proposed to describe the edge states. One picture, due to MacDonald [Q, is based upon 
a wavefunction proposed by Girvin M, which, due to particle- hole symmetry, consists of 



droplet of holes in the z/ = 1/3 Laughlin-state [ITUI embedded in a droplet of electrons in the 
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V = 1 state. This wavefunction has indeed been shown to be an excellent description of the 
exact ground-state for a system with a small number of electrons under various boundary 
conditions P,|6|,pn^ . For example, Greiter quotes an overlap of 0.9990 between the exact 



ground state and the Girvin wavefunction for a system of 8 electrons in spherical geometry. 
This z/ = 2/3 wavefunction supports two different edges 0, one at the edge of the hole 
droplet (of charge q = — e/3), and the other at edge of the u = 1 electron droplet (of charge 
q = e). On the other hand, a very different edge structure was suggested by Chang and by 
Beenakker [|12|, and elaborated on by Chklovskii et al. |jl3| in a more general context. They 



argued that for a smooth enough potential an incompressible u = 1/3 state will nucleate near 
the edges of the system, leading again to two edge branches, but this time of charge q = e/3 



each [Q. In the first scenario, where a single edge state carries a fractional charge, one 
would expect that half of the tunneling peaks will be suppressed, giving a clear signature 
of the composition of the edge states. In the second scenario, where both edges carry a 
fractional charge, all the peaks would be suppressed, resulting in a very low conductance 
signal. 

In this work we study quantitatively the nature of the ground state and the corresponding 
edge states in the v = 2/3 regime, in order to understand which of the pictures is relevant 
experimentally. Generalizing the Girvin wavefunction to incorporate the possibility of a 
V = 1/3 strip near the edge of the sample, the correlation functions in this generalized state 
are expressed exactly in terms of correlation functions calculated in the v = 1/3 Laughlin 
wavefunction. Using the mapping onto a classical one-component two-dimensional plasma 



T0| we calculate those correlation functions using classical Monte Carlo |T9| for up to 50 



electrons. The resulting z/ = 2/3 correlation functions enable us to calculate the energy of 
the state for arbitrary electron-electron interactions and confining potential. We find that as 
a function of the slope of the confining potential, the ground state makes a sharp transition 
from the Girvin-MacDonald form to the Chang-Beenakker form. This calculation suggests 
that for heterostructures where the gates are not too far from the two-dimensional electron 
gas, the suppression of half of the peaks, in the first scenario above, should be observable. 



In addition, information about the actual distance between the two edges, which is a crucial 
ingredient of recent edge state theories is obtained. The main results of this work 
have been published in a short communication [|T^. The results have also been verified by 



composite-fermion calculations M] 



The rest of the paper is organized as follows. In the following section we introduce the 
variational wavefunction and describe its construction, in Sec. Ill the correlation functions 
for this wavefunction are expressed in terms of the correlation functions for the u = 1 and 
the u = 1/3 systems, with more details in the appendix. Sec. IV describes the numerical 
evaluations of the z/ = 1/3 correlation function, while Sec.V contains the main results for 
the u = 2/3 state. Sec. VI summarizes and concludes. 

II. THE VARIATIONAL WAVEFUNCTION 

The ground state of N electrons in a radially symmetric system in the u = 1/3 fractional 



quantum Hall regime can be approximated very well by the Laughlin wavefunction |10 



where Zi denotes the complex coordinates of the i-th particle in the plane, and all lengths are 



expressed in units of the magnetic length, in = yfic/eH. = denotes a second-quantization 
representation, where a+ creates an electron in a first Landau-level state of angular momen- 
tum n, described by the single-particle wavefunction (f)n{z) = 2" exp(— |2;p/4)/-y'^27r2"n!). 
The sum is over all permutations of distinct integers which sum up to the total angular 
momentum 3N{N — l)/2, and the C^f^\„ can, in principle, be obtained by expanding the 
first product. 

The particle-hole symmetry-based wavefunction, introduced by Girvin to describe the 
u = 2/3 state consists of holes in the u = 1/3 state embedded in a first Landau level 
state of N + Nh electrons, and can be written in a second-quantized form as 

{h,--,iNf^} 



The yet undetermined number of holes, N^, must be chosen to minimize the energy. 

In order to allow for the possibility of a z/ = 1/3 state nucleating along the edge of the 
sample, we start with the Laughlin wavefunction with an inside hole of size L ||1CI|| , 

i i<J {n,---,«iv} 

(3) 

where the sum is over the same sets as in and D^'^'^^ can again, in principle, be evaluated 
by expanding the products. (Clearly D'-^'"^ = C-^^ for each set of integers.) With the above 
wavefunction describing a strip of electrons with u = 1/3 correlations along the edge of the 
sample, we write a generalized z/ = 2/3 state. 

This wavefunctionis schematically depicted in Fig.l. It depends on three integer parameters. 
Out of the N electrons, A^i are described by the Girvin wavefunction (0), with Nh holes. 
The remaining N2 = N — Ni electrons nucleate into a u = 1/3 strip along the edge of the 
sample, with minimal angular momentum L (L > A^^i + Nh). 



III. THE CORRELATION FUNCTIONS 

Our task now is to find the set of parameters Nh, Ni, and L, that minimizes the energy 
of an A^-particle state for a given confining potential and interactions. To this end we need 
to calculate the one-particle and two-particle correlation functions for this state, 

pf\z,,Z2) = Jdzs---dz^\¥''Hz„...,z^)\^ (5) 

The normalization is chosen such that pi(r) is normalized to A^, and P2{r) = 
J dzi J dz2 P2^\zi, Z2)S{\zi — Z2\^ — r^) is normalized to — 1 ||T^. The energy is easily 
obtained from these correlation functions, 



E = 2tt rdr 



V{r)p[''\r)+NUir)pi^\r) 



(6) 



where V{r) and U{r) are the potential energy and the interaction, respectively. 

In principle, of course, if one can obtain the coefficients Cj^'"' ■ i and Di^^'^^ ■ all 
correlation functions for the z/ = 2/3 state (Eq.(^) chould be readily evaluated. This, 
however, can only be achieved for a system of a very small number (A^ < 6) of particles 

= 1/3 Laughlin-like wavefunction(|3), on the other 



T8|] . Correlations functions for the 



hand, can be straightforwardly calculated for a large number of particles, using a mapping 
into a classical statistical problem ||TD|. Such a calculation, however, is not possible for 
the Girvin wavefunction (|^), or the generalized wavefunction (|), due to a resulting sign 
problem. 

The most important step in this work is expressing the correlation functions for the 
generalized z/ = 2/3 wavefunction in terms of correlation functions for the Laughlin-like 
wavefunctions for = 1/3. The main steps in this mapping are described in the appendix. 
The results are 



(2/3;N;Nf,,L,N2)/ \ 

Pi " ^ (rj 
Npt'''-^^^''-^^\zi,z,) 



Pi [r) - Pi {r) + Pi (r) 



= N,pt^\z,, z,) + N^p^i'^-^^'-Xz,, z,) + iV2pJ/^^^-^'(^,, z,) 



+ 2Re 



- 2Re 



+ 2Re 



+ Pi ('^2)Pi ('^ij-Pi (^^ijPi {r2) - Pi {r2)Pi {n, 

'3{Nh-l) Ni-1 

i=0 j=0 
L+3(Af2-l) Afi-1 

i=L j=0 
"L+3(7V2-1) 3(Af,,-l) 

Y Y <'^i>1^3'^^<^i>v3''' '^i(^i)'^*(^2)0*(2;2)</>i(2;i) 

i=L j=0 



with Ti = \zi\. The single-particle distribution function pi, is simply expressed as the sum 



of the three distribution functions of the A^^i electrons in the u = 1 state {p\ 



(l;^r-iV2)^ 



that 



of the N2 electrons in the strip of z/ = 1/3 state (p^ 



(l/3;Af2;L)^ 



and (minus) that of the 



holes in the z/ = 1/3 state, {p{ 



(l/3;Afh)^ 



. The two-particle correlation function, p2, is far more 



complicated. Nevertheless, the various terms contributing to the resulting interaction energy 
have straightforward interpretations. The first three terms describe the contribution to the 
interaction energy from interactions within the three different components. The next six 
terms describe the direct (Hartree) interactions between the three components, the Nh holes 
in the u = 1/3 state, the Ni + electrons in the u = 1 state, and the N2 electrons in the 
z/ = 1/3 strip. The last three nontrivial terms correspond to the exchange and correlation 
interactions between the different components. 



IV. THE CORRELATION FUNCTION FOR THE u = l/3 STATE 

Eq.® enables us to express the one- and two-particle correlation functions for the u = 
2/3 state in terms of quantities evaluated for the z/ = 1/3 states (Eq.|l] and Eq.^. We 
calculate the u = 1/3 correlation functions using classical Monte Carlo method |lT9| , pO| , 
based on the mapping of expectation values in the Laughlin state into statistical correlation 



functions for a two-dimensional classical plasma |]T0[, which for the wavefunction (^) takes 
the form 



exp 



y: \z.\v2 + 3 y: iog(k. - z,\') + iog(iz,n 

i i<j i 



L = corresponds to the special case Eq.(|l]). The term in the brackets can now be 
considered as the Hamiltonian for a classical system, for which Monte Carlo calculations 
can be applied. Fig. 2 depicts the one and two-particle correlation functions for 20 and 
50 particles. Similar results can be easily obtained for larger systems. As the number of 
particles increases, the edge structure in the single-particle density shifts to larger distances, 
without a significant change in the bulk density, while the two-particle correlation function 
reaches a similar maximum, but decays on a longer length scale. 

The main numerical problem in the present work is not the calculation of the correlation 
functions for large number of particles, but the increasing number of possible parameter 
sets that needs to be considered. Accordingly we limit ourselves to systems with up to 50 
electrons. 



In addition we have to calculate <ni >i/3, the average occupation of the i-th state. This 
quantity is somewhat more complicated to calculate, and so far it has been calculated with 



only a partial success . The main difficulty is that the expression of the density in terms 
of the occupation numbers, 



3(N-1) 
1=0 



(9) 



involves nonorthogonal functions, |0i(r)p. The significant overlap between them makes 
impossible the extraction of the occupations from the density profile ||2l|. In this work we 
use a novel idea to calculate the occupations using correlation functions for an — 1-state, 



r ^ 

= jdz^---dzM \{\z,-z,\'e-\^^\'l' |v^(^/^^--)(z2,...,^^)P<5(|zip-r2) 



i=2 

3(iV-l) 

= E 

where Fi can be easily obtained by expanding the Laughlin wavefunction. Comparing Eqs 



N-l) r 

J dz2--- dzN F,{Z2, ■■■,Z^) |V]/{V3;A^-) ^^^^ 



(10) 



and (|T0D shows that the occupations for the A^-electron system, < rii >i/3, can be directly 
evaluated by numerically calculating the averages of the functions Fi in the — 1-particle 
system. In Fig. 3 we show a comparison between the one-particle distribution function 
deduced from the occupations such calculated and the independently calculated distribution 
function for 20 particles. An excellent agreement is observed. In addition, the occupations 
near the edge were found to obey various exact relations obtained by directly expanding the 
Laughlin wavefunction pi] . 



V. RESULTS FOR u = 2/3 

Having obtained the correlation function for the generalized u = 2/3 state @, its energy 
can easily be evaluated for any choice of interactions and confining potential. In order to be 
as close to the experimental situation as possible, we present results for Coulomb interactions. 
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U{r) = e/er. Rather similar results have been found for other types of interactions, such as 
screened or logarithmic ones. 

As all the wavefunctions have uniform bulk density, they differ only by their structure 
near the edge, or, alternatively, by the occupations of the single-electron states near the 
Fermi energy. Accordingly, as indeed verified numerically, the only relevant ingredient of the 
confining potential is its slope at the Fermi energy. The detailed results will be presented for a 
linear confining potential, of slope S/ d, rising from zero to S over the range {rQ—d/2, rQ+d/2). 
The position of the midheight of the potential step, tq, is fixed so the filling factor is 2/3. As 
discussed below, the physically relevant parameter will be the slope of the potential, S/d. d 
is determined experimentally by the distance of the gates from the two-dimensional electron 
gas, while S is determined by the amount of voltage applied to the gate, as seen by the 
electrons in the 2d gas. For typical Ga-As samples, the gates are 120-200nm from the 2d 
gas, which corresponds to 8-12 magnetic lengths. The interaction energy e/ein is typically 
5 meV, while the boundary potential seen by the electrons is tens of meV [^. Here we will 
express all energies in units of e/ein- The calculations were done for up to 50 electrons, 
which is a typical number in an experimental quantum dot 

In Fig. 4a we plot the number of holes, Nh, which minimizes the energy for a step potential 
{d = 0), for two values of S* = 3 and S = 5. For a step potential, the ground state usually 
involves A^2 = electrons in the z/ = 1/3 strip, so it is of the Girvin type (0). The number of 
holes in the ground state is determined by the competition between the two contributions to 
the energy: the larger the number of holes, the more uniform the density, and the lower the 
interaction energy. On the other hand, the larger the number of holes, the more angular- 
momentum states are occupied, and the larger the potential energy. Thus, as the potential 
becomes softer, the number of holes may increase and a strip of electrons in the u = 1/3 
state may form near the edge. 

As can be seen from the figure, the number of holes in the ground state scales as A^/2, 
which shows the region of density different from u = 2/3 is independent of A^, namely an 
edge effect. From Fig. 4a one can obtain the physical distance between the edges for large 
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enough system, as a function of the confining potential. We find that this distance changes 
from ~ l.bin to ~ 2.5iH when 5* changes from 3 to 10. Thus, unhke the case for slowly 
varying confining potential ||13[ one cannot consider those edges as isolated from each other. 



and any effective theory should include interactions and mixing of those states . 

Since the number of holes is an integer, it will change, on average, every other time an 
electron is added to the system. This is the source of the prediction ^ that half of the 
peaks for tunneling into a z/ = 2/3 droplet will be suppressed. As the present calculation 
cannot produce the tunneling amplitudes exactly, we estimate them by their upper limit, 
the average occupation of the angular momentum state the electron tunnels to. In the inset 
of Fig. 4a we plot this occupation as a function of A^. The suppression of more than half of 
the peaks is clearly observed, with the right power-law dependence on the electron number. 
Interestingly, the calculation suggests that sometimes the ground states of and A^ + 1 
electrons differ by two holes. It remains to be seen if this is a real effect, which will result 
in a more dramatic reduction of the peak amplitude. 

In any real system the potential will rise over a finite length scale, d. We have studied 
the nature of the ground state as a function of d, and we found that for a given electron 
number A^, and potential height S, there will be a transition from the ground state being of 
the Girvin type (|^) to a state which includes electrons nucleating at the edge of the sample 
in the z/ = 1/3 state. By varying 5* it is found the transition occurs at the same ratio of 
S/d, namely at a given slope of the confining potential. For example, for 30 electrons the 
ground state evolves smoothly from = 9 and A^2 = 0, for c? = 0, to A^^ = 5 and A^2 = 
for S/d c:::^ 1, and then it changes abruptly to N^. = 15, N2 = 2 and L = Ni + + I (Fig. 
4b). Thus the two edges, of the electron droplet and the hole droplet suddenly merge, and a 
u = 1/3 strip forms, signaling a transition from the Girvin-MacDonald picture to the picture 
presented by Chang and by Beenakker and Chklovskii et al. |T^. This strip moves away 
from the edge of the electron droplet {L = Ni + Nh + 1) as S/d decreases. For example, 
for S/d = 0.6 the ground state corresponds to = 15, N2 = 2 and L — {Ni + Nh) = 20 
(Fig. 2b). For 40 electrons one can actually see two transitions. For S/c? ~ 1 the ground 
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state changes from Nh = 12 and A''2 = to A^^ = 19 and N2 = 0, namely it is still described 
by Eq.(|^), but the two edges have merged, while for 5*/^ ~ 1.4 nucleation first occurs with 
Nh = 18 and A^2 = 5. This intermediate regime where the two edges merged may suggest a 
possible description in terms of a single z/ = 2/3 edge |2^ . 

Similar transitions have been observed for other forms of confining potentials and 
electron-electron interactions. As the slope of the potential in experimental systems 
is if the order of 1.2 — 3 e^/e£^ we predict that while the experiment still being in 
the smooth side of the transition, the suppression of half of the tunneling peaks should be 
observable in quantum dots in present high mobility structures. The closer the gates to 
the 2d gas, the better the chances of seeing that effect. In addition, it is predicted that 
as a function of the voltage applied to the gates, (which changes the slope of the effective 
potential), the tunneling peak structure will change abruptly as this transition occurs. For 
high voltages half of the peaks appearing in the u = 1 regime will be suppressed in the 
u = 2/3 regime, while for lower voltages, as extensive tunneling into the = 1/3 state will 
occur, most or all of the peaks will be suppressed. 



VI. SUMMARY AND CONCLUSIONS 

In conclusion, using an exact expression for the generalized z/ = 2/3 state correlation 
functions in terms of the u = 1/3 and the u = 1 ones, we have been able to study quantita- 
tively systems of relatively large number of electrons (A^ < 50). Consequently we predict a 
transition in the nature of the ground state of the system as a function of the slope of the 
confining potential and discuss its experimental manifestation. 

In this work it was assumed that all electrons are spin polarized due to the strong 
magnetic field. There is experimental and numerical evidence that there could be a density 
and magnetic field regime, where the u = 2/3 state can involve both spin directions. The 
transition between this latter state and the spin-polarized one can be also explored using 
similar methods to the one described in this work. 
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It is hoped that this work will stimulate more experiments in this direction. 
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APPENDIX: 

Here we derive the relations between the correlation function for the z/ = 2/3 wavefunc- 
tion (Eq.^ and those for the u = 1 and u = 1/3 wavefunctions. For simplicity we detail 
the derivation for the Girvin wavefunction (Q). The generalization to the wavefunction 
is straightforward. 

We start from the definitions (^, where we expand \1/ in terms of Slater determinants 
(e.g. Eqs.(|I]) or (|^)) . To calculate pi, one can integrate over Z2, - ■ ■ , zn- The result will be 
nonzero only if — 1 of the quantum numbers of the the two Slater determinants involved 
are the same. Since the total angular momentum is fixed that means that the last quantum 
number is also identical. This leads to an equation similar to (j^), with < rii > the average 
occupations in this particular state. Since, by definition, the occupations of the Girvin 
wavefunction (|]) are trivially related to those of the full Landau level and those of the 
Laughlin wavefunction, n\^^^^ = n[^^ — n^^^^^ , the first of Eqs. (|^ immediately follows. 

The derivation of the second equation is more complicated. To calculate p2, one integrates 
over Z3, ■ ■ ■ , zn- Again, because of fixed total angular momentum, the sum of the remaining 
two quantum numbers should be the same for the two Slater determinants. Going over the 
possibihties, one finds 
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piZi,Z2) = <nin2> \(j)nAzi)\'^\(f)n2{z2)\'^ - (f)n^{Zl)K2iz2)(l)n2{Zl)(f)nAz2] 



+ <^nA^nza~ni> <P*nS^l)<P*n^{,Z2)<i)n2{.Zl)<i)nAz2). (Al) 

{ni,n2}7^{n3,n4} 

This equation is general. Next we need to relate the expectation values in the Girvin wave- 
function to the expectation values in the full Landau level and in the Laughlin wavefunction. 
A straightforward calculation gives 



< nkUm >2/3 = < nuUm >1 " < ^fc >l/3< >l - <nm >l/3< nk>l + < TlkUm >l/3 



ni "712 



n3'^4 ^2/3 



*n3 "n4 



(A2) 



where the second equation applies only for {ni,n2} 7^ {^3,?t^4}. Combining these equations 
and trivially generalizing to the wavefunction give rise to the relations (|^. 
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Figure Captions 

1. Schematic representation of the generahzed z/ = 2/3 state (^. Out of the electrons, 
A^^i are described by the Girvin wavefunction, consisting of Nh holes in the u = 1/3 state 
in the background of A^^i + Nh particles in the u = 1 state (^. The remaining N2 electrons 
nucleate into a z/ = 1/3 strip along the edge of the sample. 

2. Correlation functions for the u = 1/3 obtained from classical Monte Carlo calculations. 

3. Comparison of the density profile obtained from classical Monte Carlo calculations for 
20 particles (circles), and the density profile derived from the occupations calculated using 
Eqs.d^) and ([T0|) , and classical Monte Carlo calculations for a 19-particle system (line). 

4. (a) The number of holes that minimizes the energy of the Girvin wavefunction (0) for 
two different sizes of step potentials. The straight hues correspond to Nh oc N/2, leading to 
an N-independent edge size. The inset shows the tunneling amplitude, as estimated from 
the occupation of the relevant state, as a function of A^, on a log- log plot. The suppression 
of at least half of the peaks is evident and agrees very well with the theoretically predicted 
1/A^ dependence (b) The density profile of the ground states of = 30 electrons for 3 
different potential slopes. The existence of a u = 1 region is evident for the highest slope, 
while for the other two slopes an incompressible u = 1/3 strip is formed along the edge. 
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